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In this brief comment we show that, contrary to previous claims [Bartelt M C and Evans J W 
1993 J. Phys. A 26 2743], the asymptotic behaviour of the Kashchiev model of polynuclear growth 
is trivial in all spatial dimensions, and therefore lies outside the Kardar-Parisi-Zhang universality 
class. 
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Within the field of non-equilibrium interface growth, 
one of the central issues is the existence of dynamical 
scaling and associated universality classes |lj The 
simplest quantity which may exhibit such scaling is the 
interface width w(L,t). For a system of lateral size L, the 
scaling form for w may be written as w(t) ~ L x f (t / L z ) 
for times t larger than some microscopic time scale. In 
the limit of infinite L, one then expects the asymptotic 
time evolution of the width to follow w ~ t@ where 
f3 = xi ' z - The determination of the exponents z and 
X represents a primary objective within this field. One 
of the most popular models of interface growth is due 
to Kardar, Parisi and Zhang (KPZ) M, and through a 
concerted effort (mostly on the numerical front), there 
now exist rather precise estimates for the exponents in 
low spatial dimensions (see for example Q.) An analytic 
derivation of these exponents represents an outstanding 
theoretical challenge. 

In a recent paper || , Bartelt and Evans (BE) presented 
an analysis of the Kashchiev model ||, which is closely 
related to polynuclear growth (PNG) models @ f§ |§. 
This model has the unusual feature of being exactly solv- 
able in the sense that the average interface height h(t) 
and the width w(t) may be expressed in terms of a closed 
set of coupled integral equations. A numerical study 
of these equations was undertaken by BE. From their 
results, they concluded that the interface width of the 
Kashchiev model showed non-trivial asymptotic scaling 
with exponents which were consistent (within numerical 
error) with KPZ universality. They also raised the pos- 
sibility that an exact asymptotic analysis of this model 
may enable one to determine the upper critical dimension 
(above which (3 = 0) of the KPZ model. 

To investigate such a possibility was our motivation for 
a closer analysis of the Kashchiev model. Our results (to 
be presented below) convincingly demonstrate that the 
asymptotic behaviour of the Kashchiev model is trivial 
in the sense that the growth exponent (3(d) = 1/2 for all 
d. This model therefore lies outside the KPZ universality 
class. Our presentation shall henceforth be brief, and we 
refer the reader to ref. || for details of the formulation 
of the Kashchiev model. 

The central quantity in the Kashchiev model is the 
function 9i(t) which represents the fraction of layer i 



which is covered with deposited material at time t. The 
form of 9\(t) is known exactly, and the higher func- 
tions Oiyiit) may be generated iteratively via the integral 
equations 



, +1 (i) = / dt' {1 - exp[-(i - <r +1 ]} ^P- (1) 



df 



In terms of these layer coverages 9i, the average height 
may be expressed as 



fc(i)=X>w. 



(2) 



i=l 



and the mean square width W(t) = w(t) 2 takes the form 

oo 

W(t)=53(2i-l)0((t)-fc(t) 2 . (3) 



(These expressions differ in a minor, inessential way from 
those of BE. Note also that our unit of time is exactly 
half that used by BE.) 

The numerical analysis of BE seems to have proceeded 
by iteratively solving the set of integral equations for a 
given number of the functions 9i (t) , and then summing 
these functions in order to determine h(t) and W(t). This 
method will fail for large times as one is forced to cal- 
culate an increasingly large number of layer coverages 
to ensure numerical precision. The predictions of asymp- 
totic scaling in accord with KPZ scaling was made on the 
basis of calculating the first five layer coverages, which is 
insufficient to investigate the true asymptotic regime. 

A simple way around this problem is to sum the recur- 
sion relation (^) over all the layer coverages 9i with an 
appropriate weight such that one derives closed integral 
equations for h(t) and W(t). One then has no numeri- 
cal barriers in probing the deep asymptotic regime. The 
equation for the mean height, written in terms of the de- 
viation from linear growth A(t) = h(t) —v(d)t, takes the 
form 



1 



(4) 



l-exp(-i d+1 ) -v(d) J dt exp(-t" i+1 ). 



The d-dependent velocity is given by v{d)~ l = 
r ((d + 2)/(d + 1)), where T(z) is the gamma function 
For the mean square width, we have 



Jdf^ieM-it-tr^ 



(5) 



2h(t)- I dt' [l + 2ft(0l exp[-(t-i') d+1 ] 




1.31 t 0-1 



1 






0-1 






0.01 







5 10 15 20 25 30 



1 2 3 4 5 6 

Time t 



FIG. 1. The time derivative of the deviation of h(t) from 
linear growth, dA(t)/dt, as a function of t, for d — 1,2,3. 
Tin's data comes from numerical integration of Eq. (j^). The 
exponential decay of the envelope of oscillations is illustrated 
for the d — 3 curve. The inset shows the d-dependence of the 
decay rate X on a log-linear scale. 



The above equations are of the Volterra type, and we 
may therefore use relatively simple techniques for their 
numerical solution. More precisely, we use a uniform 
grid for the discretization, along with a trapezoidal rule 
for the integration . The discrete time step used here 
is St = 0.001, yielding results with a precision of 6 sig- 
nificant figures. We should emphasize that on reducing 
the time step further, the precision can be systematically 
improved. 

The reason for studying the deviation from linear 
growth of h(t), via the function A(t), is that one might 
expect on simple scaling grounds that A(t) ~ , In 
Fig. 1 we present our results for dA(t)/dt, for d = 1,2 
and 3. The function becomes increasingly oscillatory for 
higher dimensions. The decay of the envelope of the oscil- 
lations is exponential after some transient period (which 
grows slowly with increasing dimension.) We have deter- 
mined the decay rate X(d) for the exponential decay with 
high precision as a function of d. The inset in Fig. 1 shows 
the d-dependence of A on a log-linear scale. In Fig. 2, we 
show the evolution of the time derivative of W(t) for 
d = 1,2 and 3. Results for higher dimensions are qual- 
itatively similar. Again we see that the pre-asymptotic 
behaviour is increasingly oscillatory as one increases the 
dimension d. The asymptotic behaviour, however, is 
purely constant, implying that W(t) ~ b(d) t, and hence 
that [3 = 1/2 for all d. The full dimensional dependence 
of b(d) is plotted in Fig. 3. One may also study the 
pre-asymptotic corrections to this linear form for W(t). 
These oscillatory corrections are also found to have an 
exponentially decaying envelope, with a decay rate X(d) 
which satisfies the relation X(d) ~ 1.14(2) X(d). These 
results convincingly show that the asymptotic properties 
of the Kashchiev model are trivial; i.e. the mean height 
relaxes exponentially fast to linear growth, and the in- 
terface width evolves as the square root of time. 




FIG. 2. Tiie time derivative of the mean square fluctua- 
tions, dW(t)/dt, as a function oft, for d — 1,2,3. This data 
comes from simultaneous numerical integration of Eq. (j^) and 




FIG. 3. The d dependence of the amplitude b(d), from 
the asymptotic relation dW(t)/dt ~ 6(d). The solid line is 
the relation given in Eq. (|^), whilst the data points are from 
numerical integration of Eq. (|5j). 
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As a powerful check on the numerical work, we may 
solve the integral equation for the mean square width 
W(t) in the deep asymptotic regime by making the 
Ansatz dW(t)/dt ~ b(d), and taking h(t) ~ v(d)t : with 
v(d) defined previously. Inserting these forms into Eq. (j^) 
and taking t — > oo allows one to extract the prefactor 
b(d). Explicitly one finds 

6(d) = «(d)|r (Uf) ^) 2 -l}- (6) 

This functional form is plotted against the numerical re- 
sults in Fig. 3. Excellent agreement is obtained, indicat- 
ing both the correctness of the above Ansatz, and the 
high precision of the numerical work. 

For comparison with the above results, we mention 
that the above equations may be easily solved in the lim- 
its of d — and d — oo by use of Laplace transform meth- 
ods. One finds for d — 0, that h(t) = t and W(t) = t ex- 
actly. For d = oo one has layer-by-layer growth described 
by h(t) = n, n < t < n + 1 and W{t) = exactly. The 
infinite d results are consistent with i) the increasingly 
oscillatory behaviour found numerically for larger values 
of d (a trend towards layer- by-layer growth), and with ii) 
b(d) vanishing monotonically for large d. 

Finally we illustrate the misinterpretation of the data 
that led BE to conclude that Wit) had KPZ-type scal- 
ing. Their conclusions arose from a plot of W(t) versus 
h(t). Since their data is not in the asymptotic regime, 
such a plot is potentially misleading as h(t) still has ap- 
preciable deviations from linear growth. In Fig. 4, we 
present a linear plot of Wif) versus h(t) for d = 1, 2 and 
3 for longer times than obtained by BE (cf. Fig. 4 in ||.) 
It is seen that the relationship between W(t) and h(t) 
eventually becomes linear, in agreement with the scaling 
forms h(t) ~ v(d)t and W(t) ~ b{d)t. 

We conclude by remarking that models belonging to 
the PNG class are certainly worthy of study, as more re- 
alistic versions are known to exhibit KPZ scaling [|| 
However, any analytic treatment must be based on a less 
mean- field- like formulation than that of the Kashchiev 
model. It is of interest to determine precisely where the 
crucial mean-field assumption enters into the Kashchiev 
model, and whether it is possible to relax this assumption 
without making the model analytically intractable. 
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FIG. 4. A plot of the mean square fluctuations W(t) ver- 
sus the average height h(t), for d — 1,2, 3. 
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